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A new dark energy model called " ghost dark energy" was recently suggested to explain the 
observed accelerating expansion of the universe. This model originates from the Veneziano 
' ghost of QCD. The dark energy density is proportional to Hubble parameter, po = ctH, 

^ I where a is a constant of order Aqq^ and Aqcd ^ lOOMeV is QCD mass scale. In this paper, 

we extend the ghost dark energy model to the universe with spatial curvature in the presence 
of interaction between dark matter and dark energy. We study cosmological implications 
of this model in detail. In the absence of interaction the equation of state parameter of 
ghost dark energy is always W]j > —1 and mimics a cosmological constant in the late time, 
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^ |. while it is possible to have wd < — 1 provided the interaction is taken into account. When 
CD ■ 

■ k — 0, all previous results of ghost dark energy in flat universe are recovered. To check the 
observational consistency, we use Supernova type la (SNIa) Gold sample, shift parameter 

^ I of Cosmic Microwave Background radiation (CMB) and the Baryonic Acoustic Oscillation 

, peak from Sloan Digital Sky Survey (SDSS). The best fit values of free parameter at la 

^ '. confidence interval are: = O.Sbtom, = 0-75tg;[54 and = Omt^o^^l. Consequently 

' the total energy density of universe at present time in this model at 68% level equates to 
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I. INTRODUCTION 

The current acceleration of the cosmic expansion has been strongly confirmed by numerous and 
complementary observational data In the context of standard cosmology such an expansion 
requires the existence of an unknown dominant energy component, usually dubbed "dark energy" 
whose equation of state parameter satisfies wd < —1/3. Although we can affirm that the ultimate 
fate of the universe is determined by the feature of dark energy, the nature of dark energy as well 
as its cosmological origin is still rather uncertain, (for reviews, see e.g. [2] and references therein). 
Disclosing the nature of dark energy has been one of the most important challenges of the modern 
cosmology and theoretical physics in the past decade. A great varieties of dark energy models 
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have been proposed, to explain the acceleration of the universe expansion within the framework of 
quantum gravity, by introducing new degree of freedom or by modifying the underlying theory of 
gravity 3|-la^. 

Recently a very interesting suggestion on the origin of a dark energy is made, without intro- 
ducing new degrees of freedom beyond what are already known, with the dark energy of just the 
right magnitude to give the observed expansion [3, O]- In this proposal, it is claimed that the 
cosmological constant arises from the contribution of the ghost fields which are supposed to be 
present in the low-energy effective theory of QCD joijisl. It was argued that the Veneziano ghost, 
which is unphysical in the usual Minkowski spacetime QFT, exhibits important physical effects in 
dynamical spacetime or spacetime with non-trivial topology. The ghosts are required to exist for 
the resolution of the U{1) problem, but are completely decoupled from the physical sector [l^. The 
above claim is that the ghosts are decoupled from the physical states and make no contribution in 
the flat Minkowski space, but once they are in the curved space or time-dependent background, 
the cancelation of their contribution to the vacuum energy is off-set, leaving a small energy density 
p ~ HAqfjjj, where H is the Hubble parameter and Aqcd is the QCD mass scale of order a 
lOOMeV. With H ~ lO^^^eV, this gives the right magnitude ~ (3 x lO^^ey)'^ for the observed 
dark energy density. This numerical coincidence is remarkable and also means that this model 
gets rid of fine tuning problem [3, S] • The advantages of this new model compared to other dark 
energy models is that it is totally embedded in standard model and general relativity, one needs 
not to introduce any new parameter, new degree of freedom or to modify gravity. The dynamical 
behavior of the ghost dark energy (GDE) model in flat universe have been studied 1^ . 



In this paper we would like to extend the previous discussion on ghost dark energy [14)] to a 
universe with spatial curvature. There are enough observational evidences, at present time, for 
taking into account a small but non-negligible spatial curvature [l^. For instance, the tendency of 
preferring a closed universe appeared in a suite of CMB experiments [ij^. The improved precision 
from WMAP provides further confidence, showing that a closed universe with positively curved 
space is marginally preferred [l^]. In addition to CMB, recently the spatial geometry of the 
universe was probed by supernova measurements of the cubic correction to the luminosity distance 



18l |. where a closed universe is also marginally favored. 

Most discussions on dark energy models rely on the fact that its evolution is independent of 
other matter fields. Given the unknown nature of both dark matter and dark energy there is 
nothing in principle against their mutual interaction and it seems very special that these two 
major components in the universe are entirely independent. Indeed, this possibility has got a lot 
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of attention in the literature in recent years (see 19142 1[| and references therein) and was shown to 
be compatible with SNIa and CMB data Q . 

All above reasons, motivate us to study the interacting ghost dark energy model in a nonflat 
universe. In this paper, we would like to generalize the ghost dark energy model to the universe 
with spacial curvature in the presence of interaction between the dark matter and dark energy. 
Taking the interaction between the two different constituents of the universe into account, we study 
the evolution of the universe, from early deceleration to late time acceleration. In addition, we will 
show that such an interacting dark energy model can accommodate a transition of the dark energy 
from a normal state where wd > to wd < phantom regimes. 

This paper is organized as follows. In the next section, we review the ghost dark energy model 
in a flat universe. In section IIIH we generalize the study to the universe with spacial curvature 
in the presence of interaction between dark matter and dark energy. Observational constraints on 
the free parameters of model will be given in section IV. We summarize our results in section V. 

II. GHOST DARK ENERGY IN FLAT UNIVERSE 

Let us first review the ghost dark energy model in flat Friedmann-Robertson- Walker (FRW) 



universe where first investigated in 
to some extent from those of Ref. 



141. Although, our approach in dealing with the problem differs 



A. Noninteracting case 

For the flat FRW universe filled with dark energy and dust (dark matter), the corresponding 
Friedmann equation takes the form 

H' = ^APm + PD), (1) 

where pm and pu are, respectively, the energy densities of pressureless matter and dark energy. 



The ghost energy density is 



PD = aH, (2) 



where a is a constant of order Aq^^j-, and Aqcd is QCD mass scale. With Aqcd ~ lOOMeV and 
H ~ lO^'^^eV , A^QQYyH gives the right order of magnitude ~ (3 x 10~'^eV)^ for the observed dark 
energy density [8.]. 
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We define the dimensionless density parameters as 

Per Per 3M^H 

where the critical energy density is pcr = 3H'^Mp. Thus, the Friedmann equation can be rewritten 

as 

nm + nD = 1. (4) 

The conservation equations read 

Pm + 3Hpm = 0, (5) 
Pd + 3Hpd{1 + wd) = 0. (6) 

Taking the time derivative of relation ^ and using the Friedmann equation we find 

PD = PD^ = -^^Pd{1 + u + wd). (7) 

where u = Pm/PD is the energy density ratio. Inserting this relation in continuity equation ([6]) we 
reach 

{1 + wd){6M^H - a) = au. (8) 
Substituting ghost energy density ([2]) in Friedmann equation ([1]) we find 

3M^H = a{l + u). (9) 

Combining Eq. Q with ([8]) we reach 

= '1 + its;- 



Using the fact that 



Pm _ _ 1 — 

Pd 



(11) 



we can rewrite Eq. (|10p as 



It is easy to see that at the early time where $7/) *C 1 we have wd = —1/2, while at the late time 
where ^ 1 the ghost dark energy mimics a cosmological constant, namely u)£) = — 1. It is 
worthy to note that in wd of this model, there is no free parameter. In the left panel of figure ([1]) 



5 




we plot the evolution oi wd versus scale factor a. From this figure we see that wd of the ghost 
dark energy model cannot cross the phantom divide and the universe has a de Sitter phase at late 
time. 

We can also calculate the deceleration parameter which is defined as 

9 = -l-;^- (13) 

When the deceleration parameter is combined with the Hubble parameter and the dimension- 
less density parameters form a set of useful parameters for the description of the astrophysical 
observations. Using Eq. ^ and definition ^2/) in ([3]) we obtain 

H 3 

= --noil + u + wd)- (14) 

Substituting this relation into p^ . after using ([12]) we find 

1 3 Qp 
^^=2-2(2^ 

At the early time where rip ^ the deceleration parameter becomes q = 1/2, while at the late 
time where the dark energy dominates {^Ip 1) we have q = —1. This implies that at the early 
time the universe is in a deceleration phase while at the late time it enters an acceleration phase. 
We have plotted the behavior of q in the right panel of figure ([T]). From this figure we see that 
the transition from deceleration to acceleration take places at a ~ 0.64 or equivalently at redshift 
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z ~ 0.56. Note that 1 + z = a ^ and we have set oq = 1 for the present value of scale factor. 
Besides, taking il^, = 0.72 we obtain q ~ —0.34 for the present value of the deceleration parameter 



and using relation tlD = as well as relation (fT3|) we reach 



which is in agreement with recent observational data 

din a 

din a 

Using Eq. (fTSj) we get 



231 ] . Taking the time derivative of Eq. 



no (l + q). (16) 



= 3nD \ J" . 17 
a in a 1 — 

This is the equation governing the evolution of ghost dark energy. The dynamics of ghost dark 
energy is plotted in figure ([2]) where we have taken = 0.72 as the initial condition. This figure 
shows that at the late time the dark energy dominates, as expected. 
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FIG. 2: The evolution oiVto for ghost dark energy. Here we have taken = 0.72. 



B. Interacting case 

Next we extend the discussion to the interacting case and study the dynamics of the ghost dark 
energy. Although at this point the interaction may look purely phenomenological but different 
Lagrangians have been proposed in support of it (see [2^] and references therein). Besides, in the 
absence of a symmetry that forbids the interaction there is nothing, in principle, against it. In 
addition, given the unknown nature of both dark energy and dark matter, which are two major 
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contents of the universe, one might argue that an entirely independent behavior of dark energy is 
very special 



21 



251]. Further, the interacting dark mater-dark energy (the latter in the form of a 
quintessence scalar field and the former as fermions whose mass depends on the scalar field) has 
been investigated at one quantum loop with the result that the coupling leaves the dark energy 



2d]. Thus, 



potential stable if the former is of exponential type but it renders it unstable otherwise 
microphysics seems to allow enough room for the coupling; however, this point is not fully settled 
and should be further investigated. The difficulty lies, among other things, in that the very nature 
of both dark energy and dark matter remains unknown whence the detailed form of the coupling 
cannot be elucidated at this stage. In this case, the energy densities of dark energy and dark matter 
no longer satisfy independent conservation laws. They obey instead 




FIG. 3: Left panel corresponds to the evolution of wd for interacting ghost dark energy and different 
interacting parameter 6^ while right panel shows the evolution of the deceleration parameter for interacting 
ghost dark energy and different interacting parameter 6^. Here we took = 0.72. 



p^ + 3Hpm = Q, (18) 
Pd + 3Hpd{1 + wd) = -Q, (19) 

where Q represents the interaction term and we take it as 

Q = 3b^H{pm + Pd) = Sb^Hpoil + u), (20) 

with 6^ being a coupling constant. It is worth noting that the continuity equations imply that the 
interaction term should be a function of a quantity with units of inverse of time (a first and natural 
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choice can be the Hubble factor H) multipHed with the energy density. Therefore, the interaction 
term could be in any of the following forms: (i) Q oc HpD, (ii) Q cx Hpm, or (iii) Q oc H{pm + po)- 
Thus we can present the above three choices in one expression as Q = Tp£), where 

T = ?,b'^H for QocHpD, 

r = 3b'^Hu for QocHp,m, (21) 
T = 3b'^H{l + u) for Qo^H{pm + PD), 

It should be noted that the ideal interaction term must be motivated from the theory of quantum 
gravity. In the absence of such a theory, we rely on pure dimensional basis for choosing an inter- 
action Q. In the present work for the sake of generality, we choose the third expression for the 
interaction term. 

Inserting Eqs. ([7]) and ([20]) in Eq. (fT9]) and using (fTTj) we find 

One can easily check that in the late time where — > 1, the equation of state parameter of 
interacting ghost dark energy necessary crosses the phantom line, namely, wd = —(1 + 26^) < —1 
independent of the value of coupling constant b^. For present time with taking Q,^ = 0.72, the 
phantom crossing can be achieved provided 6^ > 0.1. This value for coupling constant is consistent 



with recent observations 



2l|. 



In the presence of interaction the deceleration parameter is obtained by substituting (j22p in 
H]) and using (fT3|) . The result is 

' ' ^ + (23) 



" 2 2(2-f]D) V 
while the evolution of dark energy follows the following equation 



din a 2 



2 — \ 17/) 



(24) 



The evolution of the cosmological parameters wd, Q and are shown in figures ([3]) and dH) for 
different interacting parameter b^. We have taken $7^ = 0.72 as the initial condition. We can also 
obtain the scale factor a as a function of t. Integrating the relation ^Id = a/{3MpH), we find 

^D-= / 17T7^dt = —,{t-to), (25) 



a 



t, 3M2 3M2 



where fi^) is given by Eq. (I24|] . The behavior of a{t) is shown in the right panel of figure 



FIG. 4: The evolution of dark energy density for interacting ghost dark energy is shown in left panel. Right 
panel corresponds to the evolution of the scale factor for interacting ghost dark energy with different . 
The rest of parameters are the same as for figure ([3]). 



III. INTERACTING GHOST DARK ENERGY IN NON-FLAT UNIVERSE 

Next we reach to the main task of the present work, namely studying the dynamic evolution of 
ghost energy density in a universe with special curvature. As we discussed in the introduction a 
closed universe is marginally favored. Taking the curvature into account, the Priedmann equation 
is written as 



where k is the curvature parameter with k = —1,0,1 corresponding to open, flat, and closed 
universes, respectively. We define the curvature density parameter as Q.^ = k/{a'^H'^), thus the 
Friedmann equation takes the form 




(26) 



1 + Ofc — Q,m + 



(27) 



Using the above equation the energy density ratio becomes 




PD 



1 + rjfc — rj^i 



(28) 



Taking the time derivative of the Friedmann equation (I26p we find 



H 



(29) 



10 



and therefore 



= PDjp = PD i^k - -^^d[^ + U + Wd_ 



(30) 



Combining this relation with continuity equation ()19p , after using ()20p and ()28p we find the equation 
of state parameter of interacting ghost dark energy in non-flat universe 

1 



The deceleration parameter is obtained as 

H 



l-^ = -l-nk + -Od[1 + u + wd] 



Substituting Eqs. ([28]) and 1^ in we obtain 



l-^ + 262fi-i(i + Ofc) 



(31) 



(32) 



(33) 




FIG. 5: The evolution of wd for interacting ghost dark energy in nonflat universe. Right panel illustrates 
the evolution of the deceleration parameter for interacting ghost dark energy in nonflat universe. Here we 
set VL^n = 0.73 and rj" = 0.28. 



In a non-flat FRW universe, the equation of motion of interacting ghost dark energy is obtained 
following the method of the previous section. The result is 



din a 



D 



1 + 



^k 



no 



^k 

3 



1 - ^ + 262^^1 (1 + Qk) 



(34) 
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FIG. 6; Left panel shows the evolution of dark energy density for interacting ghost dark energy in nonflat 
universe. Right panel corresponds to the evolution of the scale factor for interacting ghost dark energy in 
non flat universe. The value of fifc in the present time is 0.01 (closed universe). The rest of parameter are 
as in figure ([S]). 



The evolution of can be obtained by combining Eq. ^ with definition fi^t = k/[a?H'^). We 
find 

We calculated the evolution of deceleration parameter and Vto and plotted them in figures ([5]) and 
([6]), respectively. a{t) versus t in the non-flat universe for different values of coupling constant is 
shown in figure ([6]). In the limiting case Ofc = 0, Eqs. (32)- (35), restore their respective equations 
in interacting ghost dark energy model in flat universe derived in the previous section (see also 

Q). 



IV. OBSERVATIONAL CONSTRAINTS 



In this section we use the recent observational data sets for supernova type lafSN 



a) 



271, 



281], shift 



29143 ll| and Baryonic 



parameter of Cosmic Microwave Background Radiation based on WMAP-7 

Acoustic Oscillation (BAO) based on Sloan Digital Sky survey (SDSS)|32|] to put constraints on 
the free parameters of our model. To avoid the rewriting unnecessary things we refer the reader 
to some references such as j33l436l| for more details. In Table HI we summarize the list of free 



12 



TABLE I: Priors on the free parameter space. 



Parameter 


Prior 






[0.00 - 1.00] 


Top hat 




[0.00 - 1.00] 


Top hat 


Ho 




Free [37, 38] 


9 


[0.00 - 0.20] 


Top hat 



parameters of model as well as priors for using in the likelihood analysis. 
To apply the observations from SNIa we calculate the distance modulus as 



H = m-M = b\ogDL{z-A,VL%,b^) + h\og 



in the above equation 



c/Hq \ 
_ Mpc J 



dz'Ho 



+ 25, 



(36) 



H{z';no^,n%,b^)J- 



(37) 



where 



T{x) = {x, sin(a;), sinh(x)) for {^1 = 0, fl^ > 0, nl < 0) 



and H{z;fl^,Q,^,b'^) is computed numerically from Eqs. ([I]), ([2]), ([3]) and ([3 
Finally the XsNia defined by: 

il^obsizi) - fitheizi] Q°^, n%, 62)]2 



(38) 



Usually, beside using the peak locations of the CMB power spectrum, one can use the so-called 



shift parameter TZ, as 



221 



7^= aAH 



■DLizdec,^L^%,b^ 



(39) 



(1 + Zdec) 

here Zdec is the redshift of the last scattering surface [40]. Subsequently the Xcmb ^e written 



as 



2 



(40) 



For the last observational constraint, we rely on the large-scale correlation function measured 
from the sample of SDSS including a clear peak at 100 Mpc/i~^ 32]. A dimensionless and Hq 
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TABLE II: The best fit values for tlie free parameters from fitting witli SNIa from new Gold sample, 
SNIa+CMB, SNIa+CMB+BAO experiments at one and two a confidence level. 



Observation 






62 


SNIa 


0.95_o.56 

0.95_o.8o 


0.99-0.25 
0.99-0.44 


0.16-0.11 
0.16-0.16 


SNIa+CMB 


n "30+0.22 
"•"-•^-0.16 

U.O»_g 3Q 


"• '"J-0.17 

n 70+0.08 
"• ' "^-0.26 


nq+o-io 

0.09-0.09 


SNIa+CMB+BAO 


"•>53-0.07 


n 7^+0-01 
n 7S+0-04 

"• '"J-0.08 


o.o8;^:°^3 
0.081^:0^5 



independent parameter for constraining the cosmological models has been proposed in literatures 
3^ as follows: 



A 



where z^^^g — 0.35 



1/3 



32l |. So the Xbao expressed as: 



(41) 



(42) 



BAO 



Figure ([7|) represent the marginalized likelihood function for model free parameters. In addition 
joint contour plot for parameters have been illustrated in figures ([8]) and Q 

The best values and the confidence interval for free parameter at la and 2a have been reported 
in Table dnl). 



V. CONCLUSION 



It is a general belief that our universe is currently undergoing a phase of accelerated expansion 
likely driven by dark energy. Unfortunately, until now, the nature and the origin of such dark 
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FIG. 7: Marginalized likelihood functions of model free parameters. The solid, dash and dashdot lines 
correspond to fitting the model with SNIa data new gold sample, SNIa+CMB and SNIa+CMB+BAO, 
respectively. The horizontal solid and dashed lines represent the bounds with la and 2a level of confidence, 
respectively. 



energy is still the source of much debate and we don't know what might be the best candidate 
for dark energy to explain the accelerated expansion. Thus, various models of dark energy have 
been proposed, to explain the accelerated expansion by introducing new degree of freedom or by 
modifying the standard model of cosmology. In this regard, a so called "ghost dark energy" was 
recently proposed [3, 0] which originates from the Veneziano ghost of QCD. The QCD ghost has 
no contribution to the vacuum energy density in Minkowski spacetime, but in curved spacetime it 
gives rise to a small vacuum energy density [8]. The dark energy density is proportional to Hubble 
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FIG. 9: Joint likelihood function of free parameters, left panel corresponds to {b'^,il'^) while right panel 
shows for {b^,V,^). Here SNIa+CMB+BAO observations used to confine the values of free parameters. 



parameter, pD = aH, where a is a constant of order Aq(-.j-) and Aqcd is QCD mass scale. With 
Aqcd ~ WOMeV and H ~ lO^^Sg^ ^ 

gives the right order of mag nitude ^ (3 x lO^^eV)^ 
for the observed dark energy density [8]. The advantages of this new proposal compared to the 
previous dark energy models is that it totally embedded in standard model so that one needs not 



4ll |. Besides, some experimental data has 



16 

to introduce any new parameter, new degree of freedom or to modify general relativity {l^ . 

In this paper, we generalized the ghost dark energy model, in the presence of interaction between 
dark energy and dark matter, to the universe with spatial curvature. Although it is believed that 
our universe is spatially flat, a contribution to the Friedmann equation from spatial curvature is 
still possible if the number of e-folding is not very large 
implied that our universe is not a perfectly flat universe and recent papers have favored the universe 
with spatial curvature [3]. With the interaction between the two different dark components of the 
universe, we studied the evolution of the universe, from early deceleration to late time acceleration. 
We found that in the absence of interaction the equation of state parameter of ghost dark energy 
is always larger than —1 and mimics a cosmological constant in the late time. We also found that 
the transition from deceleration to acceleration take places at a ~ 0.64 or equivalently at redshift 
z ~ 0.56. We observed that, in the presence of interaction, the equation of state parameter can 
cross —1 at the present time provided the interacting parameter satisfy > 0.1. 

To check the observational consistency of interacting Ghost Dark Energy model, we used Super- 
nova type la (SNIa), CMB shift parameter and Baryonic Acoustic Oscillation (BAO). Our results 
demonstrated that the best values of free parameters when we combine all observational data are: 
17^ = 0.35+[j:[!^, n% = 0.75toM and 6^ = 0.08lg;|j| at la confidence interval. Our analysis shows 
that at la level of confidence the value of so-called interacting parameter does not cross zero. Also 
the total value of energy density of universe at present time is fi^Qt = + $72) = LIO^J^qq^ at 68% 
level. 

Finally, we would like to mention that if there is any kind of ghost field which gives rise to 
an energy density p (x H, its cosmological implications is exactly similar to the present work 
independent of its origin. Although the existence of a well-motivated physical model where this 



energy density is obtained is of course a valic 



be found in the previous works such as and we have not repeated them here. In this work 



starting point, however, the details of this study can 



our main tas' 
proposed in 



i was to study the cosmological implications of this new ghost dark energy model 



i 



without referring to its origin. In particular we generalized the study to the 



universe with spatial curvature and put some observational constraints on the model parameters. 
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